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Abstract. It is widely accepted that the low-energy physics in edge-sharing cuprate materials 
has one-dimensional (ID) character. The relevant model to study such systems is believed to 
^ ■ be the ID extended Heisenberg model with ferromagnetic nearest-neighbor (NN) interaction 

and antiferromagnetic next-nearest-neighbor one. Thus far, however, theoretical studies of such 
materials have been confined to the case of isotropic interactions. In the present work, we 
compare the spin susceptibility of the ID extended Heisenberg model with anisotropy in the 
NN channel, obtained by means of the Transfer Matrix Renormalization Group method, with 
£3 . that of the edge-sharing chain cuprate Li2Cu02. 
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1. Introduction 

In the last decades, an increasing attention has been paid to low-dimensional materials showing 
frustrated magnetism. It has been found, by using several experimental techniques, that the 
low-energy physics in such edge-sharing cuprate materials has one-dimensional (ID) character 
and develops in the chains of CUO4 plaquettes. A spin 1/2 residing on a Cu + ion interacts 
with its nearest neighbor via an oxygen ion in such a way that the Cu-O-Cu bond forms an 
CT\ I obtuse angle and, therefore, the nearest-neighbor (NN) exchange interaction among Cu + ions 

appears to be small and ferromagnetic. The interaction between next-nearest neighboring Cu + 
ions develops on a Cu-O-O-Cu path via the exchange mechanism and is antiferromagnetic. On 
the basis of these considerations, it is widely accepted that the effective model for edge-sharing 
cuprate materials is a ID extended Heisenberg model (EHM) with ferromagnetic NN interaction 
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and antiferromagnetic next-nearest-neighbor (NNN) one: 
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In order to estimate the parameter values for the effective EHM model, electronic-structure 
calculations and experimental data from susceptibility and specific heat are usually combined 
with theoretical curves obtained by means of either Exact Diagonalization or Transfer Matrix 
Renormalization Group (TMRG). 

A comparison between temperature-dependent quantities, such as susceptibility, specific heat 
and magnetization, and theoretical predictions for model ([TJ has already been done for e.g. 
Pb[CuS04(OH)2] as in Ref. [HE] or for Rb2Cu2Mo30i2 as in Ref. [3]. Such comparisons showed 



a good agreement at high temperatures while, at low temperatures, the necessity to include 
inter-chain coupling, which could bring the system into the 3D regime at T —* 0, or additional 
Hamiltonian terms, like Dzyaloshinskii-Moriya, has been argued. 

Up to now, all attempts to compare theory with experiments have been done for the isotropic 
model ([1]) only. However, the anisotropy of the interactions in such systems has already been 
emphasized as for LiCuV04 [4] and for Li2Cu02 [5]. From the theoretical point of view, 
anisotropy lowers the symmetry of the system and could radically change the properties of 
the spin model already at zero temperature [UJ [SJ E] . 

In the present work, we compare the magnetic susceptibility of L12CUO2 compound, taken 
from Ref. [TO], with that of the ID extended anisotropic Heisenberg model 

H = -J Z J2 S?S? +1 + J ± J2(SfSf +1 + SfS? +1 ) +J'J2 SA+2, (2) 

i i i 

calculated within TMRG and try to infer the anisotropy ratio J±/J z by fitting the experimental 
curves. For the sake of simplicity we assume anisotropy only in the NN channel, since in Li2Cu02 
NN coupling is dominant. 

2. Method 

A detailed description of the TMRG method can be found in [11] . In successive diagonalization- 
decimation of the transfer-matrix, we have used the Arnoldi algorithm with implicit restart |12j 
followed by the bi-orthogonalization of the left and right eigen-spaces of the non-symmetric 
density matrix. In doing the decimation, we retain at most 100 lowest eigen-states of the 
density matrix. We have verified that the truncation error in our calculations was less than 
10~ 4 . Since the model ([2]) contains also the NNN term and the TMRG is designed for systems 
with only the NN one, we adopt the mapping of the Hamiltonian ([2]) proposed in Ref. [13]. 

3. Results 

On the basis of the zero-temperature phase diagram obtained for the model in Refs.JSJE], 
we expect Li2Cu02 to be in a particularly interesting region of the phase diagram In fact, it 
could fall either into a massive phase or into a massless one, depending on the precise values 
of Jj_ and J'. Our strategy is as follows: since the anisotropy estimated from the experimental 
data is expected to be small [HE], we first choose J' in order to fit best the high-temperature 
behavior, then introduce a small deviation of Jj_ from J z and choose the value of J± which fits 
best the experimental susceptibility curves at intermediate temperatures. In Ref. [10j . we have 
found three sets of data for the susceptibility of Li2Cu02, corresponding to the three possible 
orientations of the magnetic field along the principal axes, as shown in Fig. [IJ being the 6-axis 
along the chain. 

There is a long and controversial history related to the estimation of coupling constants 
J and a = J' / J in Li2Cu02, for which the isotropic model (UJ) is usually assumed and hence 
J = J z = J±. It was originally believed, after Ref. [TO], that in Li2Cu02 J = WOK and a = 0.62. 
Lately, first from quantum chemistry calculations [2] and then from exact diagonalization of 
pci-Hubbard model, representing finite chains of Cu02 plaquettes [15j . different estimates came 
out: J = 142K, a = 0.15 in the first case and J = 146K, a = 0.23 in the second one. 
Finally, from Density Functional Theory calculations [16] the values J = 215K and a = 0.31 
were claimed. We have explored the range of a z = J'/J z £ [0.15,0.62] and NN anisotropy 
a_i_ = J_l/J z € [0.6,1.3]. The effect of the anisotropy on the susceptibility at fixed a z is to 
decrease the height of the intermediate-temperature peak upon increasing a±. Since for a± < 1 
this peak is overestimated, (not shown in Fig.l because its scale is too different) the range 
a±_ < 1 is ruled out and we end up with a± > 1. We have found that the only set of a z and 
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Figure 1. First three panels: theoretical fits for magnetic susceptibility of model ([2]) to the 
experimental data for Li2Cu02 at a z = J±/J z = 0.33, different values of a± = J±/J z and for 
different orientations of magnetic field. The values of ^-factors, extracted from Curie's law for 
different field orientations are: g a = 2.15, = 1.92, g c = 1.85. J z = 215K and J' = 0.33J Z 
were taken. Last panel: zoomed view in the vicinity of the transition temperature T c of the 
experimental data and theoretical curves for all three orientations of magnetic field. 



a± compatible with the experimental data in the widest possible range of temperatures is the 
following: a z = 0.33 and a± = 1.25. We have also determined the values of g-factor by requiring 
that, in the high-temperature limit, the theoretical susceptibility asymptotically tends to the 
corresponding experimental one. This gave us the following values: g a = 2.15, = 1.92 and 
g c = 1-85. 

Below T c = 9.1K, it was found experimentally that in Li2Cu02 a long-range order, 
ferromagnetic along the chain (6— axis) and antiferromagnetic between the chains (a— axis) 
is stabilized [T7]. Hence, even a small but finite inter-chain antiferromagnetic coupling will 
dramatically reduce the magnetic response at low temperatures. That is why our theoretical 
susceptibility (see Fig.l last panel) overestimates the magnetic response of Li2Cu02 at T < T c . 



4. Conclusions 

In the present manuscript, by using TMRG, we have investigated the possibility to improve the 
theoretical description of magnetic susceptibility in the edge-sharing cuprate material Li2Cu02 



upon introduction of a small anisotropy in the NN channel. The main effect of the NN anisotropy 
at fixed J z appears to be a change in the height of the susceptibility maximum at intermediate 
temperatures. Namely, the increase of anisotropy a± suppresses the maximum, while the 
decrease of a± has the opposite effect. Within the whole range of possible values of a z , known 
from the literature, we found that a z = 0.33 together with the anisotropy a± = 1.25 best 
describe the experimental susceptibility data in terms of the model ([2]). We believe that taking 
into account the inter-chain coupling along a— axis would further improve our description for 
the corresponding susceptibility at T < T c and such work is currently in progress. 

Acknowledgments 

The numerical calculations reported in the present article were done in part on CINECA CLX 
cluster (project No. 954) and on the CASPUR cluster (project No. 315/09). 

References 

[1] Kamieniarz G, Bielinski M, Szukowski G, Szymczak R, Dyeyev S and Renard J P 2002 Comp. Phys. Comm. 
147 716 - 719 

[2] Baran M, Jedrzejczak A, Szymczak H, Maltsev V, Kamieniarz G, Szukowski G, Loison C, Ormeci A, Drechsler 

S L and Rosner H 2006 Phys. Stat. Sol. C 3 220-224 
[3] Lu H T, Wang Y J, Qin S and Xiang T 2006 Phys. Rev. B 74 134425 (pages 9) 

[4] Krug von Nidda H A, Svistov L E, Eremin M V, Eremina R M, Loidl A, Kataev V, Validov A, Prokofiev A 

and Afimus W 2002 Phys. Rev. B 65 13445 
[5] Ohta H, Yamauchi N, Nanba T, Motokawa M, Kawamata S and Okuda K 1993 J. Phys. Soc. Jpn. 62 785-792 
[6] Plekhanov E, Avella A and Mancini F 2008 Acta Phys. Pol. A 113 429 
[7] Plekhanov E, Avella A and Mancini F 2008 J. Optoelectron. Adv. M. 10 1675 
[8] Avella A, Plekhanov E and Mancini F 2008 Europhys. J. B 66 295 
[9] Plekhanov E, Avella A and Mancini F 2008 preprint cond/mat-0811.2973 
[10] Mizuno Y, Tohyama T, Maekawa S, Osafune T, Motoyama N, Eisaki H and Uchida S 1998 Phys. Rev. B 57 

5326-5335 

[11] Xiang T and Wang X Density-Matrix Renormalization: A New Numerical Method in Physics ed Peschel I, 

Wang X, Kaulke M and Hallberg K (New York: Springer 1999) p 149 
[12] see http : / /www . caam . rice . edu/sof tware/ARPACK 
[13] Maeshima N and Okunishi K 2000 Phys. Rev. B 62 934-939 

[14] de Graaf C, de P R Moreira I, Bias F, Iglesias O and Labarta A 2002 Phys. Rev. B 66 014448 
[15] Malek J, Drechsler S L, Nitzsche U, Rosner H and Eschrig H 2008 Phys. Rev. B 78 060508 
[16] Drechsler S L, Malek J, Nishimoto S, Nitzsche U, Kuzian R, Eschrig H and Rosner H 2009 J. Phys. Conf. Ser. 
145 012060 

[17] Sapina F, Rodriguez-Carvajal J, Sanchis M, Ibianez R, Beltran A and Beltran D 1990 Sol. State Comm. 74 
779 



